Parity Effects in Spin Decoherence 



in 
o 
o 

< 



> 
o 

(N 

O 

Q-i! 
^— > ■ 
C ■ 

=5 

cr 



X 



A. Melikidze, 1 V. V. Dobrovitski, 2 H. A. De Raedt, 3 M. I. Katsnelson, 4 and B. N. Harmon 2 

1 Kavli Institute for Theoretical Physics, University of California, Santa Barbara CA 93106, USA 
2 Ames Laboratory, Iowa State University, Ames I A 50011, USA 
3 Applied Physics - Computational Physics, Materials Science Centre, 
University of Groningen, Nijenborgh 4, NL-9747 AG Groningen, The Netherlands 
' Uppsala University, Department of Physics, SE-751 21, Uppsala, Sweden 
(Dated: April 14, 2004) 

We demonstrate that decoherence of many-spin systems can drastically differ from decoherence 
of single-spin systems. The difference originates at the most basic level, being determined by parity 
of the central system, i.e. by whether the system comprises even or odd number of spin-1/2 entities. 
Therefore, it is very likely that similar distinction between the central spin systems of even and odd 
parity is important in many other situations. Our consideration clarifies the physical origin of the 
unusual two-step decoherence found previously in the two-spin systems. 

PACS numbers: 03.65.Yz, 75.10.Jm, 76.60.Es, 03.65.Ta 



INTRODUCTION 

Reduced dynamics of a small quantum system coupled 
to a bigger environment has recently become the subject 
of particularly active investigation. In fields like quan- 
tum optics [lj and quantum computation, Q, there is a 
naturally defined distinct "central" system (i.e. an atom 
or a qubit) which interacts with its environment, and 
whose dynamics is of primary importance. Similar sit- 
uations are often encountered in the condensed matter 
physics, e.g., when considering a heavy particle tunnel- 
ing in a crystal, tunneling centers in glasses Kondo 
systems etc. This problem is also of importance when 
a naturally defined central system is absent, such as in 
a recently developed promising approach to the theory 
of strongly correlated systems, the dynamical mean- field 
theory (for review see |j|). In this approach, the sys- 
tem of interacting particles in a crystal is replaced by an 
"effective impurity" in a self-consistently defined thermo- 
stat. 

So far, quantum evolution of a single two-level sys- 
tem (or, equivalently, a single spin-1/2 entity) interact- 
ing with a bath of bosons 5] or spins 0, Q has been 
studied in much detail. In contrast, the central sys- 
tems comprising several strongly interacting spins 1/2 
have not been that extensively investigated. A general 
analysis of the two-spin central system interacting with 
a bath of bosons has been presented in Ref. , but more 
detailed considerations are lacking. Several interesting 
cases of a two-spin system coupled to a spin bath have 
been considered in Refs. 0,0, and it has been demon- 
strated that behavior of many-spin central systems can 
be very different from a single-spin case. Consideration 
of many-spin central systems is of particular importance 
for possible implementation of quantum computations 
which use several strongly coupled two-level systems for 
encoding of a single qubit 0, 0]. This representa- 
tion allows using the "decoherence-free subspaces" and 



error-correcting schemes developed for multi-spin qubits 

BUI. 

In this work, based on an exactly solvable but realis- 
tic model, we show explicitly that decoherence of a two- 
spin-1/2 system can be qualitatively different from deco- 
herence of a single spin 1/2. We demonstrate that this 
difference originates at the most basic level, and is de- 
termined primarily by parity of a central system, i.e. by 
whether the central system comprises even or odd num- 
ber of spin-1/2 entities. 

It is known that the parity of the spin system is the 
cause of the drastically different behavior in the tunnel- 
ing of magnetization in a wide class of spin systems such 
as magnetic nanoparticles and molecular magnets where 
the tunneling is due to magnetic anisotropy or magnetic 
field 0,0]. In this paper, we explore a different effect, 
in which the parity of the central system determines the 
long-time dynamics of the decoherence process. We em- 
phasize that in the system considered here the quantum 
oscillations are caused by the isotropic exchange interac- 
tion and are independent of the symmetry of the crystal 
field and external magnetic field; thus the short-time os- 
cillations do not depend on the parity. 

Although there are many possible central systems cou- 
pled to various kinds of spin baths, the generic differ- 
ences between the many-spin and the single-spin central 
systems can be understood based on simple models. An 
instructive model of a many-spin central system inter- 
acting with a spin bath, has been recently analyzed by 
Dobrovitski et al. 3] This model is aimed to describe 
(at least, qualitatively) main features of such central sys- 
tems as magnetic molecules, quantum dots or impurity 
spins which experience decoherence from the nuclear spin 
bath. In these systems, the dominant interaction with 
the nuclear spins can be approximated by the isotropic 
Heisenberg interaction, since anisotropic interactions are 
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often small. The model is defined by the Hamiltonian: 

N 

H = C 2 + 2CY / JkSk, (1) 

fe=i 

which describes the central system composed of two 
spins: C — Ci + <?2, Ci = C2 = 1/2, which is coupled 
by Heisenberg exchange interaction to N environmental 
spins Sfc = 1/2, k = 1 . . . N. Note that the environmen- 
tal spins don't have their own dynamics. This may be 
viewed as a limit case where the dynamics of the central 
system is much faster than that of the environment. 

A special feature of this model, which makes it differ- 
ent from the "central spin" models considered by Garg, 
or Prokof 'ev and Stamp, ||J, b| is the fact that in our 
treatment the central system is not reduced to the dou- 
blet of lowest states. This features is crucial to the results 
discussed below. 

One is interested in the time-evolution of the ini- 
tial system-plus-environment state which is taken in the 
form: 

N 

K)=iT) cl u) C2 ni^- ( 2 ) 
*i=i 

The initial states of the environmental spins \i) Sk are as- 
sumed random and uncorrelated. The initial state of the 
system is a superposition of the singlet and triplet states 
of the two central spins: 

lt)c 1 |l) C2 = ^(|l,0)c + |0,0)c), (3) 

where we have introduced notation \C, C z )c for the cen- 
tral spin. One considers the problem of the decay of this 
coherent singlet-triplet superposition in the central sys- 
tem due to its interaction and subsequent entanglement 
with the environmental spins. In particular, one is in- 
terested in the time-dependence of the expectation value 
of the z-component of the first spin (of (i)), where of is 
the Pauli matrix acting on the state of c\. In the ab- 
sence of the coupling to the environment this quantity 
exhibits periodic oscillations between +1 and —1 caused 
by the first term in Eq. ijljl: coupling to the environment 
is expected to damp these oscillations. 

In the work reported in Ref. [|| a numerical investiga- 
tion of this problem was performed. Among many sur- 
prising features in the behavior of the above system, it 
was observed that after an initial fast decay of the oscil- 
lations of (af(t)) the amplitude showed a saturation at 
the value of 1/3. Subsequently, the oscillations demon- 
strate a much slower decay, which is consistent with the 
\jt conjecture, and which leads to a complete suppres- 
sion of oscillations. The main motivation of this paper 
was to understand the cause of the saturation and the 
subsequent slow decay. 



While the model Eq. is hard to treat analytically, 
we simplified it by setting all J^'s equal. This allowed 
us to solve the model exactly. The solution turned out 
to reproduce quantitatively several key features of the 
numerical results reported in 0. In fact, it reproduced 
the fast initial decay of the amplitude of oscillations and 
its subsequent saturation at 1/3. It also offers a way to 
qualitatively understand the cause of the long-time tail. 
Most importantly, it answers the question: why is the 
decay of oscillations in our model much slower compared 
to a more conventional exponential decay of oscillations 
in, say, the spin-boson models. [f| The cause is the integer 
value of total spin of the central system. 

This work shows that integer spins, in contrast to half- 
integer spins, may, under suitable circumstances, exhibit 
quantum oscillations over much longer times. From the 
perspective of the theory of quantum phase transitions, 
this work also offers a simple example of emergent power- 
law correlations usually associated with criticality. 

J fc = J MODEL 

To make analytical progress we consider a simplified 
model where we take all coupling constants Jfc — J to be 
equal while preserving random uncorrelated initial states 
of the environmental spins. The Hamiltonian takes the 
form 

H = C 2 + 2JCS 

= (1-J)C 2 + J{C + S) 2 ~JS 2 , (4) 

which describes the coupling of the central spin C — 
c\ + c 2 to the total spin of the environment S = J2sk- 
We are interested in the expectation value of the z- 
component of c%: (of (i)), where of is the Pauli matrix 
acting on the state of c\. Note, that the assumed ini- 
tial condition Eq. J2J corresponds to the superposition 
of states with different S. The Hamiltonian Eq. con- 
serves S 2 , therefore the matrix element (of(t)) can be 
decomposed as 

K(i))=£<S|<7^)|S)P(S), (5) 
s 

where P(S) is the weight of the state with the total spin 
S given the random uncorrelated initial states of We 
thus are led to the problem of first calculating P(S). 

Before proceeding with the actual calculation an im- 
portant comment is in order. Since Eq. (JjJJ looks like an 
average over all possible initial orientations of the envi- 
ronmental spins one might interpret the above quantity 
(of(t)) as an ensemble-averaged expectation value. Quite 
importantly, in the case where the number of environ- 
mental spins is large the actual weight of the state with 
total spin of the environment S tends to the ensemble- 
averaged quantity P(S). Therefore, in this limit Eq. (J5J 
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describes well the evolution of the central system in a 
single realization of the experiment. 

In the basis where Sk are good quantum numbers the 
initial density matrix is, by assumption, a 2 N x 2 N ma- 
trix: 



-N 



I. 



(6) 



where / is a unit matrix. Let us make a unitary trans- 
formation to the basis spanned by the the eigenstates of 
S 2 . There are N different values that 5 12 can take. To 
preserve the dimensionality of the Hilbert space we con- 
clude that some (in fact almost all) of these latter states 
are degenerate. A unitary transformation will leave the 
initial density matrix unchanged. This means that 



P(S) 



-N 



G(S)(2S+1), 



(J) 



where G(S) is the degeneracy of the state with total spin 
S (with S z fixed). To calculate G(S) we change vari- 
ables and introduce g(k) = G(N/2 — k). The state with 
the maximum total spin S = N/2 is unique and is the 
state where all s^'s point up (we choose S z — S), there- 
fore g(0) = 1. Next, a state with S = S z = N/2 - 1 
should be a superposition of the states with N — 1 spins 
up and one spin down. There are C± — N such states 
(C*^ = N\/M\{N-M)\ is the binomial coefficient). How- 
ever, among such states there are g(0) — 1 states with 
S = N/2 and S z = N/2 - 1 which have to be excluded. 
Generalizing to arbitrary k we get: 



fc-i 



g(k) = c, 



JV 



k ~^2,9{i) 

i=0 

\ ,N-2k + l 



'h 

fit 
L 'k 



N 



1 



(8) 



We thus have the result for the weight of the state with 
spin S: 



P(S) 



N/2-S 



(2S+1) 2 

N/2 + 5+1 



8S 2 /N 
V2nD 



-S 2 /2D 



D = N/4, (9) 



where we have used a well-known approximation for the 
binomial distribution described by the first two factors 
above. One can easily check that J Q P(S)dS = 1, i.e. 
the approximations we made preserve the normalization 
of the probability. 

We have thus reduced the problem to finding the time 
evolution of the initial state: 



If) 
I*) 



e- im \i) 



V2 



(|l,0)c + |0,0)c)|S,^)s, 



(10) 
(11) 



estimating the spin polarization (f\<rf\f), and averaging 
the result with respect to S z (trivial) and S (according 
to Eq. jnj). There are two circumstances that greatly 
simplify the calculation. First, the Hamiltonian acting 
on the singlet state |0,0)c gives zero, therefore the evo- 
lution of the second term in Eq. I|ll|l is trivial. Second, 
the symmetry of the Hamiltonian with respect to c\ and 
<?2 implies that given the above initial condition we have 
(/|<rf |/) = — (/lofl/)- Thus, we can calculate the expec- 
tation value of ci = {a\ — c|)/2 instead. For this op- 
erator we have: ot|l,0) c = |0,0)c, ot|0,0) c = |l,0)c, 
crf_ j 1, ±l)c = 0. Taking all this into account, we see that 



<JK|/) = Re(i|e-^), 
\t) = \l,0) c \S,S')s. 



(12) 
(13) 



From Eq. (J2J it is clear that the above matrix element 
can be easily calculated after going to the basis with well 
defined total spin L = C + S. The necessary Clebsch- 
Gordan decomposition (in the limit S ^> 1 of interest to 
us) is: 



'l- (S z /S) 2 



+ -g\S,S*) L , 



(\S + 1,S Z ) L -\S-1,S Z ) L ) 

(14) 



where we have introduce the notation |L,L z )i. In this 
basis we easily calculate using Eq. 10} : 



(fWllf) = cos2(l- J)t 



' cz \ * fa 

< [1 - ( — I lcos2JS't 



S 



S 



(15) 



Finally, we have to average this result over S z and S. The 
first average is done trivially using the fact that (in the 
same limit S > 1) ({S z ) 2 ) = S 2 /3. The second average 
is calculated using Eq. (O which leads to a Gaussian 
integral. The result is: 



K) = A(t)cOs2(l - J)t, 

A(t) = l + l(l-NJH 2 )e 



(16) 
(17) 



It should be stressed that this result is exact in the limit 
N ^> 1 (S ^> 1). We see that an initial exponential 
decay of the amplitude of the oscillations is followed by 
a transient and an eventual leveling at A(t) = 1/3. 

To check the above results we have performed a di- 
rect numerical solution of the Schrodinger equation corre- 
sponding of the system with a Hamiltonian H = JqG 2 + 
2CJ2k=i JkSk, Jk = J, which can be reduced to Eq. 
by rescaling J — > 1, J — * J/ Jo, t — * tJ . Exact diagonal- 
ization was used to find the time evolution. An example 
of the results is shown in Fig. 1. It shows the expectation 
value of erf as a function of time. The parameters are: the 
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Time 

FIG. 1: Numerical simulation of 13 spins with Jo = 8, Jk — 
J — 0.128. The figure shows the expectation value of af as a 
function of time. 




5 10 15 20 

Time 

FIG. 2: Analytical result for af (t) with the same parameters 
as those used in numerical simulations. 

number of spins N = 13, J = 8, J k = J = 0.128. This 
can be compared with the analytic result for the same 
quantity which is given (after rescaling) by Eq. (|l(jfl and 
is shown in Fig. 2. The numerical and analytical results 
show excellent agreement. 

Absence of the decay of the amplitude of oscillations 
at long times is quite an unexpected result. Therefore it 
is worth explaining it in some more detail. 



DISCUSSION OF THE RESULTS: SIMPLE 
PICTURE 

One trivial situation where the oscillation of the central 
spin does not decay is that of no interaction between the 
central spin and the set of environmental spins. In the 



presence of such an interaction, however, one may still 
ask what are the conditions under which this interaction 
is ineffective in damping the oscillations. A natural sug- 
gestion is to try to find a state |^) of the combined system 
in which (*|i7 int |*) = 0. Since in our case H int = 2JCS, 
classically such a state |\E') would correspond to vectors 
C and S being orthogonal. The condition CS = defines 
a plane in 3D space, therefore one could argue that for 
the case of random initial orientation of S the probability 
of being in the state with H lnt = is zero. Remarkably, 
the quantum nature of spins proves the result to be quite 
different. 

The correct way of treating H mt is, of course, to rewrite 
it in the following form: 

flint = J{C + S) 2 - JS 2 - JC 2 . (18) 

Adding spin C = 1 with spin S results in possible values 
of total spin C + S being: S — 1, S and 5+1. It is the 
second case in which the first two terms in the Eq. Ijl8(l 
cancel each other. The remaining last term does not 
depend on S and, therefore, does not suppress the oscil- 
lation amplitude when the averaging over S is performed 
and only shifts the oscillation frequency of the central 
spin (this effect is reflected in Eq. I jlfijl ). The condition 
"C + S = S" is, thus, the closest analog of the classical 
condition CS = 0. But, unlike in the classical case, a sim- 
ple Clebsch-Gordan algebra (see previous section) shows 
that the probability of being in the subspace "C+5 = S" 
is actually finite and is equal to 1/3. 

One can easily see now that this effect can only oc- 
cur if the central system has integer spin. Indeed, the 
condition "C + S — S" can never be satisfied if C is half- 
integer. These considerations allow us to formulate the 
main result of the paper: Based on a particular model 
of a central spin interacting with randomly-oriented en- 
vironmental spins we have been able to show that the 
decay of the oscillations of the central spin is essentially 
different for integer central spins: the decay is no longer 
exponential, instead the amplitude of the oscillations sat- 
urates at a constant value. 

Moreover, the results presented in this work make clear 
the physical origin of the unusual two-step decoherence 
found in Ref. |2j, where the generic model Eq. Q has 
been considered with all Jk being different. The first step 
of decoherence, associated with the initial decay of oscil- 
lations to the value of 1/3, has been described in Ref. 
using a mean-field-like treatment of the spin bath, by 
replacing the interaction part of the Hamiltonian with 
a random classical static field having Gaussian distribu- 
tion. However, such a treatment fails to describe the 
second step of decoherence, i.e. the long-time slow decay 
of oscillations. As the results above demonstrate, the 
representation of a bath as a static random field corre- 
sponds to the case of all Jk being equal to J. This stems 
from the fact that the total spin of the bath S 2 commutes 
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with the Hamiltonian @, so that the bath dynamics in 
the case Jk = J is trivial, and can be removed completely 
by a transformation into the rotating coordinate system. 
Then, in the rotating coordinate system the effect of the 
bath on the central spins is equivalent to the action of 
a random static field. Therefore, the initial decoherence 
is similar to the "adiabatic decoherence" by a static spin 
bath, considered e.g. in Ref. |l7j . 

Correspondingly, the second step of the decoherence 
process, i.e. the long-time slow decay of quantum oscilla- 
tions, can be caused only by an internal evolution of the 
bath. For all Jk being different, S 2 does not commute 
with the interaction part of the Hamiltonian (@J , and, as 
a result, the system-bath coupling induces a non-trivial 
dynamics inside the bath. It is not surprising that the 
spin bath possessing a complex dynamics can not be rep- 
resented as a random static magnetic field acting on the 
system. Understanding this "minimally non-adiabatic" 
decoherence regime represents a challenge for future in- 
vestigations. pj| 

Summarizing, in this work we have demonstrated that 
decoherence of many-spin systems can drastically differ 
from decoherence of single-spin systems. This difference 
originates at the most basic level, and is determined by 
parity of the central system, i.e. whether the system com- 
prises even or odd number of spin-1/2 entities. Therefore, 
it is very likely that similar distinction between the cen- 
tral spin systems of even and odd parity is important in 
many other situations. Moreover, our consideration clar- 
ifies the origin of the unusual two-step decoherence found 
numerically in Ref. Q . The exactly solvable model allows 
clear demonstration that the initial step of decoherence 
(associated with the saturation of oscillations at the value 
of 1/3) is caused by "adiabatic decoherence" by a static 
spin bath, while the subsequent long-time slow decay is 
induced by a non-trivial internal dynamics of the spin 
bath. The model is applicable to the qualitative analy- 
sis of a range of experimental systems such as magnetic 
molecules and shallow impurity spins in semiconductors, 
which experience decoherence from the nuclear spin bath. 
In these cases, the dominant interaction with the nuclear 
spins is well approximated by the isotropic Heisenberg 
interaction (anisotropic interactions are often small). 
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